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Abstract 

We show that the Moyal ★-product on the algebra of fields induces an effective 
lattice structure on vortex dynamics which can be explicitly constructed using 
recent asymptotic results. 
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1. Introduction 

The study of the behavior of classical fields defined as functions of noncom- 
mutative spatial variables has received a great deal of attention in the last few 
years (see e.g. and Q for a review). In particular the study 

of coherent structures in the form of noncommutative solitons or noncommu- 
tative vortices has shown that the noncommutative version of ip'^-type models 
in two spatial dimensions with polynomial nonlinearities sustain non-collapsing 
soliton plateau and uncharged static vortex solutions, unlike what occurs in 
the commutative case. It has also been shown that the confining mechanism 
is provided by the ★-product which acts like a projection. Using this fact it 
has been further shown that large vortices can be produced by taking static 
suitable combinations of projections. In this process vortex-like structures with 
no angular dependence appear, and which have quantized radii R ~ ^yn for n 



'Corresponding author 
Email addresses: timOmym . iimas .unam.mx (A.A. Minzoni), 
roman. juarezSnucleares .unam.mx (L.R. Juarez), mrosenOnucleares. unam.mx (M. 

Rosenbaum) 



Preprint submitted to Elsevier 



February 7, 2009 



integer. It is also known QjQl that spatial noncommutativity in 2-spheres in 
R"^ induces a lattice structure in the radial variable, where the quantized radii 
scale as i? ~ n. Moreover this lattice can confine radial solutions. 



On the other hand, recent work in non-linear optics (see e.c 



B. 0. 



12l | and [13[ for a review and background) has shown that optical lattices can 
trap coherent structures in the form of plateaus and vortices, with and without 
charge. The dynamics of these structures has been described successfully using 
modern asymptotics in terms of the Peierls-Nabarro (P-N) potential 1^ , 14 1 
produced by the interaction between the coherent structure and the physical 
lattice. 

The purpose of this letter is to show by means of an asymptotic analysis how 
the ★-product combined with the vortex structure induces a lattice in the spatial 
dimensions via a P-N like potential in the radial variable which, in turn, confines 
the vortex itself avoiding the collapse. We show, making use of a coherent state 
vortex-type solution with varying parameters, how the average Lagrangian [isl l 
for the noncommutative Nonlinear Schrodinger equation is equivalent to the 
known average Lagrangian for a vortex on a discrete lattice. The difference 
between the two Lagrangians being the actual form of the P-N potential. For 
the lattice the P-N potential is periodic with the lattice period, while on the 
noncommutative case considered here the P-N potential depends on R^, which is 
the square of the vortex radius. This dependence gives a scaling R ^ ^Jn which is 
expected from the exact result described for the (p'*-type models. We also exhibit 
how for vortices in two and three spatial dimensions, whose charge is of the same 
order of their large radius, the ★-product induces a P-N potential which results 
in an equispaced radial lattice. This shows how the splitting of in terms of 
fuzzy spheres placed on an equispaced radial lattice, as proposed in QjBj arises 
naturally from ★-product confining a vortex with large charge, large radius and 
small width. This results show how the dynamics of the coherent structures 
sustained by the noncommutativity can be described asymptotically in terms 
of the classical dynamics of coherent structures in lattices. Furthermore, the 
present analysis shows that the qualitative behavior for large vortices does not 
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depend on the details of the noncommutativc model chosen. 



2. Formulation 

We consider the Nonlinear Schrodinger equation where the nonlinearity is 
given by the Moyal ★-product between the complex functions u(x, t), x is the 
2-dimensional spatial position and t is the time. 

The equation is given by the usual local Lagrangian [is! ]. suitably modified 
by the nonlocal product in the form: 



C = J J d xdt[i{utu — utu) + [^ul — (1) 

where the Fourier transform of the ^-product of two complex fields u and v in 
the Rieffel formula for the Moyal ^-product [16] is given by: 

5^(k) = ^-i^ / u{k - p)w(p)e*(''-P)^Pdp. (2) 
(27r)^ J 
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Here {t(k, i) ~ fe ^^'^u{K.,t)d:!C and aAb = §(ai,a2) 

with Q fixing the square of the length given by the noncommutativity of the 
spatial variables. 

The equation associated with the Lagrangian ^ is the nonlocal Nonlinear 
Schrodinger equation: 

iut = Au + 2u-k u-k u. (3) 

This is the usual Schrodinger equation where the real potential [/(x) is now 
dependent on the solution u-k u. 

In previous studies 1] of the Sine-Gordon equation, static (with no angular 
momentum) exact vortex- type solutions were found in the limit of — > cxd. 

Analogous solutions of ^ can be obtained in the form 

u(x,t) = e'-*C(x), (4) 

where ?; is a real valued function of the position. Substitution into ^ readily 
gives 

-aC- AC-fC*C*C (5) 
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In the large noncommutativity limit ^ becomes the same as the equation con- 
sidered in [l| . Specially interesting solutions to ([5]) in this range of values of the 
noncommutativity parameter discussed in the above cited work are given (using 
scaled variables) by: 

(^„(r) = 2(-l)"e-'-'L„(2r2), (6) 

with r'^ — + and cr = 40 and where L„ are the Laguerre polynomials. It 
is known that sums of ^Pn{f) are also solutions. For example, taking 

n 

Pn{r) = ^ (fijir) and Wn = Vn(^) + V'n-iir). 
j=o 

Thus, using the well established fact that the Laguerre polynomials have caustics 
when ^ n [13], we obtain that the plateau P,i{r) is practically flat up to 
r ^ ^Jn and confined to that region. On the other hand, Wn has a peak at 
r ^ ^fn. Such a behavior suggests trapping by an annular lattice with radii 
^ ^/n, for large n. 

Moreover, it is known both from analytical as well as numerical cal- 

culations that the usual Nonlinear Schrodinger equation on the usual square 
lattice supports stable vortices due to the trapping of the vortex by the P-N po- 
tential, which prevents their collapse as it occurs in the limit of the continuum. 
We will show below, using the asymptotic analysis developed in [l^, that the 
★-product in equation ([3]) indeed generates, via the field, a P-N like potential 
responsible for the trapping of the vortex. We will further show how this P-N 
potential generates lattices with radii growing as 0i or as n, depending on the 
coherent state used for the averaging of the Lagrangian. 



3. Asymptotic solutions 

On the basis of the above considerations, we shall derive next an asymptotic 
solution to 13]) in the form of a vortex with angular momentum (and charge one) 
given by the local behavior re*^ at the origin, and an envelope suggested by the 
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exact solution with no angular momentum. Thus, we take as the coherent state 
trial function to average the Lagrangian ([T|) the expression 

u{r,e,t) = a(i)re-(^)%^(^+'^)e'('^-«(*))^W, (7) 

where the amplitude a, the width co and the phase a are functions of time. The 
radius R represents the location of the peak of the vortex and the velocity V is 
the radial velocity of the maximum of the vortex. We will also consider charged 
vortices with charge m where the angular dependence is r™e™^ close to r = 0. 

Let us begin by studying the charge one vortex by substituting the trial 
function ([7]) into the Lagrangian ([1]) and performing the spatial integration we 
obtain an averaged Lagrangian for the parameters a, lu, R, V and a. 

This Lagrangian is then varied and the modulation equations obtained for 
the parameters which give the approximate evolution of the vortex. The aver- 
aged Lagrangian has the form 

C = Co + C,, (8) 

where Cq arises from the local terms in ([1]) and is the contribution of the 
★-product. We have, using the results in J4], that 

= -2c?ujR^(j - 2a^R^uj{VR - — ) . (9) 

To calculate we begin by transforming the potential term with the 
product in ([T]) into the Fourier space. We have, using the Parseval relation, 
that 



-ku-ku-kud^ — I \u-ku\^dx^ / \u -k u\'^ (k) dli. (10) 



The calculation of the last integrand is performed by making use of and 
the Fourier transform of coherent state ([7]) which is given by: 



OC f27T 



u(p) = a(t)e*"(*M / e'f"°^(^-^)e*^re-(^) e'('-«'^r(i(^dr, (11) 



/o Jo 

where p = p(cos(6'), sin(6')) and (p is the polar angle of x. 



5 



In the highly noncommutativc hmit the width uj is small. Hence the inte- 
grand in r is peaked at r = i?. Making the change of variables r = i? + lj^ we 
obtain, to leading order in w, 



/ e'P^''°''^e"^e''"^^''°''^e-^\"^^^rd^dr. (12) 

-R Jo 

Since we are interested in large vortices, we use the stationary phase approx- 
imation [l8^ on the angular integral. There are two stationary phase points at 
if — 0, If ^ n which give an oscillatory contribution. The radial integral is then 
calculated extending the lower limit to minus infinity to obtain: 

(27r)3/2pi 
In the same way we obtain 



Hp) - {p" + j) (13) 



«(P) - J0^''- ("«+ l) (14) 

Consequently, using (|13p and (fT4|) together with ([2]) allows us to arrive at an 
approximation for the integrand in (llOp in the form 



(.-„)(.) . .^.v f i ' ^^^^ 

X sin [pR + ^ sin {\k - p|i? + ^) e'''^e'''^pdpdeid02, 

where 9i and 62 correspond to the angular coordinates of p and k — p, respec- 
tively. Again (|15p will be evaluated approximately in the strongly noncommu- 
tativc limit using the stationary phase method in the angular integral. To this 
end recall that both k and p are large since the vortex is narrow because co is 
small. Using p = — we obtain 



kAp= -fcgsin6l, |p-k| = W-^ -2-^/ccos6'-Hfc2 (16) 

Since in the strongly noncommutativc limit k is also large we know that the 
points of the stationary phase are for 9 — tt/2 and 9 — 37r/2. Using again the 
same type of calculation as in the derivation of equation \12^ we obtain: 
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IjJ I \ V / / 



X sin ( —q ) sin ( R\\ + fc^ ) ^ (t'^dq. 

This integral is again evaluated asymptotically for larger i? using the method 
of stationary phase. We need to observe that for small and to leading order 
of this width the stationary phase point is g = . We then have 



X sin sin2(fci?)e-^'='«i+^)'+5(i-^)'). 



(18) 



Finally, integrating this last expression over k we obtain that the Lagrangian 
term is given by: 

(19) 

and, evaluating once more with the method of stationary phase for large R we 
obtain: 



(20) 



where the function F is the analogue of the Peierls-Nabarro potential for the 
lattice generated by the noncommutative self interaction of the field, and takes 
the form: 



Hence the final average Lagrangian is: 



1/2 

£ =2a'^ujR^& + 2a^ujR + 2a^(jR^{VR - — ) 
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This Lagrangian is, except for the dependance in the potential, the same 
average Lagrangian obtained in [cite] for large Nonlinear Schrodinger commu- 
tative vortices on a discrete lattice. However in the present case the lattice is 
generated by the ik-product and it manifests itself in the dependance of the 
Peierls-Nabarro potential. 



4. Vortex solutions and radial stability 

The approximate dynamics of the vortex is obtained from the variational 
equations of P^ . These are 



Sa : ^{2a^ujR) = 0, 

5a : AaujR^— + — — + daF = 0, 



dt 


3a; 


da 






3t^2 



5oj : 2a^R'—^—^ + d^F = 0, (23) 
dt Suj"' 

5V : R-V - dvF = 0, 

6R : ^2a^ujR^V + OrC = 0. 
dt 

The last expression above is the equation of motion for the peak of the vortex, 
analogous to a particle in the Peierels-Nabarro potential. 

The dynamics of the solutions to simplifies for large vortices in lattices, 
with low kinetic energy, i.e. for V <^ 1. In this case, since u is small, the 
dominating terms for the 5a and equations can be readily solved to give a 
steady vortex amplitude width relation in the form 

a = 32yf. (24) 

For large vortices with Rlo still small and with small kinetic energy V \ 
the dynamics of the peak is described by the simple equation 



8 




(25) 



Equation shows that the vortex moves in the lattice which was generated 
by the ★-product. The fixed points are the possible equilibrium positions and 
are given by 



The odd values of n give stable vortices, while the even values of n give insta- 
bility. If a vortex starts at an unstable value it will shrink radiating until it is 
trapped at the lower minimum of the potential. It is to be noted that the scaling 
of the radius is the same as the one obtained by the exact trapped solutions 
of (p^-type noncommutative models. This result has a simple interpretation in 
terms of the Landau cells used in [tIi. In fact the P-N potential generated by 
the T^-product induces an annular lattice, where the n-th annulus has an area 
An = 2TrRn{Rn+i — Rn) and where _R„ is given by (|25p . As n — > cx) so that 
An ~ TT^Q, which gives the constraint area of the Landau cell. We thus can say 
that the lattice induced by the ★-product is a lattice of Landau cells. The same 
calculation when performed for a charge m vortex, with R of the same order as 
m, gives an equispaced lattice. In fact the asymptotic evolution of the integrals 
in Eq.(I12]) replaces the term sin^(fci?) of HH) by sin^(fc + l)R. This results in 
a P-N potential of the form (PT|) where the term (cos + ^) is replaced by 
cos(_R + f ). this induces an equispaced lattice. The area of the Landau cells is 
again constant, as can be readily verified. Finally, when we go to R"^ and use 
in the averaging of the Lagrangian coherent structures with no angular depen- 
dence, we obtain shell lattices with i?„ ~ y^; while when including an angular 
dependence given by the spherical function Yim{9, ip), induces - for large values 
of I and as a consequence of the form of the coherent state - equispaced lattices, 
thus recovering in a natural way the noncommutative structure proposed in Q . 



Rn = V2mTQ. 



(26) 
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5. Discussion and Conclusions 

We have shown that the effect of noncommutativity of the spatial variables, 
when averaged on the appropriate coherent vortex or plateau-like states, induces 
an effective spatial lattice of Landau cells whose distribution and sizes depend on 
the coherent states in question. This shows that the effect of noncommutativity 
on coherent structures whose local width is comparable to the spatial scale Q of 
the 7k-product behave as classical structures on a physical lattice and allows us 
to calculate the lattice and the corresponding dynamics of the noncommutative 
coherent structure. It is to be remarked that the lattice structures in three space 
dimensions have been constructed using a group on the sphere and taking the 
Casimir values as the possible values of the radii inducing a lattice. We have 
shown that in the appropriate coherent states the :*r-product induces the same 
lattice in the radial variable. 

We also observe that unlike physical lattices which are not translation in- 
variant, the lattices induced by the ^-product are translation invariant. Because 
of this reason coherent states can move with uniform velocity. We thus see that 
the effect of the ^-product has no classical analogy since it is capable of forming 
a lattice to support the coherent structure without loosing the uniform transla- 
tional motion. 

We end by remarking that in the weak noncommutative limit, which is the 
opposite limit to the one considered here, equation ([T]) resembles a high order 
non-linear system of equations which incorporates Raman scattering effects [l^ . 
However at the present time the existence of an optical analogue of equation ([1]) 
in the strongly noncommutative limit is not known. 
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